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ABSTRACT 

\q • We present and analyse light curves of four transits of the Southern hemisphere extra- 

ly-) \ solar planetary system WASP-4, obtained with a telescope defocussed so the radius 

. of each point spread function was 17" (44 pixels). This approach minimises both ran- 

fT) ■ dom and systematic errors, allowing us to achieve scatters of between 0.60 and 0.88 

£. | , mmag per observation over complete transit events. The light curves are augmented 

■ by published observations and analysed using the JKTEBOP code. The results of this 
' process are combined with theoretical stellar model predictions to derive the physical 
, properties of the WASP-4 system. We find that the mass and radius of the planet 

are M b = 1.289 ±S:§io ±8:800 M Jup and R h = 1.371 +g;™§ +°;°ou Rju P , respectively 
.J^ I (statistical and systematic uncertainties). These quantities give a surface gravity and 

density of g h = 17.03 ±g;|| ms" 2 and p h = 0.500 ±8:8i? ±0.008 l°Jup» and fit the trends 
$-h ' for short-period extrasolar planets to have relatively high masses and surface gravities. 

WASP-4 is now one of the best-quantified transiting extrasolar planetary systems, and 
significant further progress requires improvements to our understanding of the physical 
properties of low-mass stars. 

Key words: stars: planetary systems — stars: individual: WASP-4 — stars: binaries: 
eclipsing 
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1 INTRODUCTION 

The d iscovery of the first extrasolar planet l|Mavor fc Quelo3 
1995) heralded the arrival of a new field of astronomical re- 
search: the investigation of how planets form and evolve. To 
date over 350 extrasolar planets have been discovered, most 
through the radial velocity motion of their parent stars. The 
resulting sample is statistically valuable, but only limited 
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information can be extracted for individual objects. The ex- 
istence of transiting extrasolar planets (TEPs) provides the 
solution to this problem, as for these objects it is possible to 
determine their masses, radii and temperatures, and hence 
surface gravities, densities and to some extent chemical com- 
positions. At present about fifty TEPs are known, but the 
discovery rate is increasing fast. 

Our ability to determine the physical properties of a 
TEP is very dep endent on the qua l ity of the available ob- 
servational data. ISouthworthl J2OO8, 2009) presented an ex- 
haustive study of fourteen TEPs, to provide an homoge- 
neously measured set of physical properties which can form 
the basis of statistical studies of these objects. This work 
highlighted the fact that measurements of the properties of 
TEPs depend critically on the quality of their transit light 
curves. The limiting factor for such observations is that the 
transits are very shallow (only a few percent of the star's 
light is lost) so it is only possible to measure transit shapes 
well if observat ional errors are very low . 

In Paper I (jSouthworth et alj2009h we investigated tele- 
scope defocussing as a way to minimise the random and 
particularly the systematic errors which afflict astronomical 
time-series photometry. This technique has a long history 
but has only very recently been regularly applied to charge- 
coupled device (CCD) observations (see references in Pa- 
perl). The most important advantage is that flat-fielding 
errors decrease according to the square-root of the number 
of pixels in a point spread function (PSF), so the effects 
of flat-fielding can be lowered by orders of magnitude com- 
pared to focussed observations. Changes in atmospheric see- 
ing and telescope pointing affect photometry through flat- 
fielding errors, so are similarly diminished. Another benefit 
is that longer exposure times are possible without saturat- 
ing the CCD, which can thus be read out less often. The 
main downside is that more background light is collected, 
but for many situations this is unimportant. The signal-to- 
noise calculations presented in Paper I show that the opti- 
mum defocus, in the sense of achieving the best signal to 
noise per unit time, is surprisingly large and can lead to 
exposure times of many hundreds of seconds even during 
periods of high sky brightness. The limiting factor is actu- 
ally the need to adequately sample the temporal variation 
of transit events. 

In this work we present defocussed photometric ob- 
servations of four transits of the Southern TEP WASP- 
4, which was discov ered by the Super WASP collaboration 
(Wils on et al.l 120081). This object has since been studied 
bv lGillon et al.l (|2009h . who used t he Very Lar g e Tele scope 
(VLT) to observe one transit, and IWirm et al l (|2009h . who 
measured two transits with the Magellan Baade telescope. 
These two studies both used defocussing to improve their 
photometric precision, and the larger telescope apertures 
(8.2 and 6.5 m) also allowed a relatively high observing 
cadence. Our observations used a 1.5 m teles cope, have a 
photo metr ic precision compa rable to those of iGillon et al.l 
(2009) and|W mn et al.l(|2009f ). and cover more transit events 
but at a lower sampling rate. In this work we will analyse 
th ese three sets of tran sit observations, plus those obtained 
by I Wilson et al.l (|2008f > using the 1.2 m Euler telescope. Our 
meth ods are identical to those used bv ISouthworthl (2008, 
2009), so are fully homogeneous with the physical properties 
of fourteen TEPs determined in these works. 




Figure 1. Surface plot of the PSF of WASP-4 in an image taken 
at random from the observing sequence on the night of 2008 Oc- 
tober 1st. The x and y axes are in pixels. The lowest and highest 
counts are 587 and 25823 electrons, respectively, and the z axis 
is on a linear scale. 

2 OBSERVATIONS AND DATA REDUCTION 

We observed four transits of WASP-4 using the 1.54 m 
Danish Telescope at ESO La Silla and the DFOSC focal- 
reducing imager. This setup yielded a full field of view of 
13.7'xl3.7' and a plate scale of 0.39" pixel -1 . The CCD 
was windowed down to 1200x1000 pixels to decrease the 
readout time from approximately 90 s to 30 s. An observing 
log is given in Table[T] All observations were done through 
the Cousins R filter, using exposure times of 120 s. The 
amount of defocussing was adjusted until the peak counts 
per pixel from WASP-4 were roughly 25 000 above the sky 
background, resulting in a doughnut-shaped PSF with a di- 
ameter of about 44 pixels (17"). The pointing of the tele- 
scope was maintained using autoguiding, and we did not 
change the amount of defocussing during an observing se- 
quence. An example PSF is shown in Fig.[T] 

Several images were also taken with the telescope prop- 
erly focussed, and were used to verify that there were no 
faint stars within the defocussed PSF of WASP-4 which 
might dilute the transit depth. We found that the closest 
detectable star is much fainter and at a distance of 31" (78 
pixels), so the edge of its PSF is separated from that of 
WASP-4 by 34 pixels in our defocussed images. The closest 
star of similar to or greater brightness than WASP-4 lies 80" 
away. We conclude that no stars interfere with the PSF of 
WASP-4. 

Data reduction was performed in the same way as in 
Paper I, so we only give a summary here. We used a cus- 
tom pipeline written in the i DiJfl programm ing language and 
using the daophot package (|Stetsonll 19871 ) to perform aper- 
ture photometry with the APER routine (distributed as part 



1 Information on the 1.54 m Danish Tele- 
scope and DFOSC can be found at 
http : //www. eso . org/sci/f acilities/lasilla/telescopes/dlp5/ 
^ The acronym IDL stands for Interactive Data Language and 
is a trademark of ITT Visual Information Solutions. For further 
details see http://www.ittvis.com/ProductServices/IDL.aspx 
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Table 1. Log of the observations presented in this work. N obs is the number of observations and 'Moon' is the fractional illumination of 

the Moon at the midpoint of the transit. 

Date Start time (UT) End time (UT) N ^ s Exposure time (s) Filter Airmass Moon Scatter (mmag) 

2008 08 19 04:47 10:24 134 120.0 R c 1.11 -» 1.03 -* 1.54 0.936 0.84 

2008 08 23 04:30 10:19 126 120.0 R c 1.10 -> 1.03 -» 1.59 0.580 0.88 

2008 09 23 01:15 05:04 88 120.0 R c 1.27 -> 1.03 -» 1.04 0.402 0.75 

2008 10 01 00:57 05:27 103 120.0 R c 1.23 -» 1.03 -> 1.10 0.035 0.60 
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Figure 2. Final light curves of WASP-4 from our four nights of 
observations. The error bars have been scaled to give x 2 „ = 1.0 
for each night, and are mostly smaller than the symbol sizes. 

of the ASTROLl4f| library). Apertures were placed by eye 
then fixed throughout each observing sequence. The aper- 
ture which gave the most precise photometry was of radius 
26 pixels, with a sky annuhis of 33-50 pixels; the shape of 
the light curve is not sensitive to the aperture sizes. 

Five comparison stars were measured on each im- 
age, and subsequently checked for short-period variability. 
Many differential-magnitude light curves were calculated for 

3 The ASTROLIB subroutine library is distributed by NASA. For 
further details see http://idlastro.gsfc.nasa.gov/ 



Table 2. Excerpts of the light curve of WASP-4. The full dataset 
will be made available in the electronic version of this paper and 
at the CDS. 



HJD 


Differential magnitude 


Uncertainty 


2454697.700127 


0.00040 


0.00086 


2454697.702257 


0.00055 


0.00085 


2454701.688470 


0.00052 


0.00085 


2454701.690241 


0.00018 


0.00085 


2454732.553043 


0.01181 


0.00079 


2454732.554814 


0.01570 


0.00079 


2454740.540454 


-0.00117 


0.00049 


2454740.542225 


-0.00061 


0.00049 



Table 3. Literature times of minimum light of WASP-4 and their 
residuals versus t he ephemeris deriv ed in this work. 
References: (D lGillon et alj <2009l) ; (2) IWinn etafl <2009l) ; (3) 
This work. 



Time of minimum (HJD) Cycle O — C value Reference 

number (HJD) 



2453963.1086I ;oo2i -549.0 0.00021 1 

2454364.5757±q o°33 -249.0 -0.00214 1 

2454368.592661°;^^ -246.0 0.00012 1 

2454371. 26738tg;°ggg7 -244.0 -0.00162 1 

2454396.69548t°;°g°26 - 225.0 0.00008 1 

2454697.797489 ± 0.000055 0.0 0.00000 2 

2454748.650490 ± 0.000072 38.0 0.00021 2 

2454697.79748 ± 0.00010 0.0 -0.00001 3 

2454701.81218 ± 0.00013 3.0 -0.00000 3 

2454732.59118 ± 0.00013 26.0 -0.00033 3 

2454740.620821 ± 0.000061 32.0 -0.00008 3 



WASP-4 and between comparison stars. Slow variations in 
brightness were observed for almost all of these, and are at- 
tributable to atmospheric effects. We therefore defined time 
intervals outside transit events and fitted straight lines to 
the magnitudes to rectify the light curve shape. Simultane- 
ously with this procedure, the comparison stars were com- 
bined into a weighted ensemble which produced the final 
light curve with the least scatter in the intervals outside 
transit. We found that including flat- field corrections pro- 
vided only a slight improvement in every case, as the tele- 
scope pointing was good, and that accounting for the CCD 
bias had a negligible effect. We have applied bias and flat- 
field corrections to generate our final light curves, which are 
shown in Fig. [2] and reproduced in Table(2] The scatter in 
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the final light curves varies from 0.60 to 0.88 mmag, and the 
most precise data were obtained in dark time when the sky 
background was low. 



3 LIGHT CURVE ANALYSIS 
3.1 Period determination 

The first step in defining the physical properties of the 
WASP-4 system is an accurate orbital ephemeris. We have 
obtained individual times of mid-transit by analysing the 
data from our four nights separately, using the jktebop 
code and Monte Carlo simulations (see below). These have 
been supplemente d by the two transi t timings given by 
IWinn et all (|2009l) and the five listed bv lGillon etahl (|2009 ) 
found from their 'prayer-bead' analysis. 

We have assigned cycle numbers to the times of min- 
imum light, taking as reference epoch the midpoint of our 
first transit (which by coinci dence is the same t ransit event 
as the first one observed by I Winn et al.l [2009h . and fitted 
a straight line to the transit timings. We find an orbital 
ephemeris of 

T = HJD 2 454 697.797488(31) + 1.33823150(61) x E 

where E is the number of orbital cycles after the reference 
epoch and quantities in parentheses denote the uncertainty 
in the final digit of the preceding number. The reduced x 2 
of the fit is Xu = 1-35, which suggests the possibility that 
the orbital period is not constant. The transit timings are 
collected in Table0 along with the residuals of the straight- 
line fit, which plotted in Fig. [3] 

In the course of measuring the times of minimum light 
given by our data, we found that the absolute sizes of the ob- 
servational error bars delivered by our pipeline (which orig- 
inate from the aper algorithm) are not reliable. We have 
rescaled the e rror bars by mult i plying them by \/ xS t o 



1.00 



give Xv =1-0 l|Bruntt et alJl2006l : [Southworth et alj[2007al ). 
treating each night individually. 



3.2 Light curve modelling 

We have analysed the light curves using an identica l ap- 
proach to that discussed in detail bv lSouthworthl (|2008l ). We 
therefore retain homogeneity with the fourteen TEPs anal- 
ysed in that work, as well as with our study of WASP-5 (Pa- 
per I). The_ligW_cyrrveswere_modelled using the JKTEBOfjfl 
code liSouthworth et all l2004al lbh . which is based on the 
ebop pr ogram originally developed for eclipsing binary star 
syste ms (|Popper fe Etzellll98ll ; lEtzellll98ll ; rNelson fc Davisl 
1 19721 ). This simulates the two components using biaxial 
spheroids, so allows departures from sphericity. The param- 
eter which governs the shapes of the components is the mass 
ratio; we have adopted the value 0.0013 and find that even 
large changes have a negligible effect on the solution. 

The main results from the light curve analysis are the 
best-fitting orbital inclination, i, and fractional radii of the 
star and planet, va = ^ and rb = — S where R\ and Rb 
are the absolute radii of the components and a is the orbital 



4 JKTEBOP is written in FORTRAN77 and the source code is avail- 
able at http://www.astro.keele.ac.uk/~jkt/ 
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Figure 4. Phased light curves of WASP-4, compared to the best 
fits found using JKTEBOP and the quadratic LD law with both 
LD coefficients included as fitted parameters. Top is from the 
Danish telescope (t his work); second is f rom the Euler tele scope 
iWilson et al.l2008l); third from the VLT llGillon et 3.1.1120091) ; last 
is from Magellan Twinn et al.ll2009h . The residuals of the fits arc 
plotted at the bottom of the figure, offset from zero. 



semimajor axis, ta and are actually parameterised as 
their sum and ratio, (rA+J"t>) and k — —, because these 
two quantities are only weakly correlated for a wide variety 
of transiting planet and other eclipsing binary systems. 

We included the limb dark ening (LD) of the star us- 
ing five different LD laws (see ISouthworthl |200S| for their 
definition) with the coefficients of these laws either fixed at 
theoretically predicted values, included as fitted parameters, 
or perturbed around theoretical values. The theoretical co- 
efficients w ere obtained by bilin e ar inter p olation in the tab- 
ulations of IVan Hamme] lll993h . IClaretl (|2000l . l2004al i and 
IClaret fc Hauschildtl l|2003l ). 

The uncertainties of the fitted parameters were as- 
sesse d using 1000 Monte Carlo simu lations for each solu- 
tion (|Southworth et al.ll2004c] . l2005bh . We have also tested 
for correl ated ('red') noise u sing the residual-permutation 
method (| Jenkins et alj [2002). finding a small but signifi- 
cant amount which is probably of instrumental origin. How- 
ever, it is possible that it is caused by spots on the sur- 
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Figure 3. Plot of the residuals of the times of mid-transit of WASP-4 versus a linear ephemeris. Some error bars are smaller than the 
symbol sizes. 



face of the parent star, based on the shapes of our transit 
light curves at HJDs 2454701.79 and 2454732.61 (Fig.0. 
The other two tra nsits, one of which was also observed by 
W inn et al. I (|2009T I, are seemingly unaffected. 

At the request of the referee we have calculated the 
parameter j3, a measure of the systematic noise contribu tion 
(|Winn et al.ll2007l : iGillon et aLlbood : IWinn et al.ll2008l ), for 
comparison with other works. We evaluated /3 values for 
individual transits and for groups of between two and ten 
datapoints. The average values of /3 over these groups are 
1.03, 1.40, 1.41 and 1.14 for our four observing sequences. We 
emphasise that we do not use these numbers in our analysis, 
but account for systematics using the residual permutation 
approach instead. 

The final light curve fits consist of sets of best- fitting pa- 
rameter values with Monte Carlo uncertainties, for each LD 
law and for different combinations of fixed/perturbed/fitted 
LD coefficients. In each case the linear coefficient is repre- 
sented by ua and the nonlinear coefficient by da- The solu- 
tions for all LD coefficients fixed are given in Table[4] and 
for all LD coefficients fitted in TableJS] The latter solutions 
are clearly better, indicating that the theoretical LD coef- 
ficients are not a perfect match to the data. We therefore 
adopt these solutions, with uncertainties which have been 
increased to account for the correlated noise discussed in 
the previous paragraph. A fit to the data is shown in Fig. [4] 
for which the rms of the residuals is 0.85 mmag. 

In order to maximise the accuracy o f the final results, w e 
have also analysed the data obtained bv lGillon et al.l |2009), 
I Winn et al. (|2009 | ). an d the Euler Telescope data presented 
by I Wilson et al.l l|2008h . The full sets of solutions of all four 
light curves are given in the Appendbfl (Tables |AH lAH lAljl 
and IA4|) . In each case we adopted the results for both LD 
coefficients fitted, and performed Monte Carlo and residual 
permutation analyses to obtain the final errors. The best fits 
are shown in Fig.|4] and the final parameters for each light 
curve are collected in Table[B] 

Also in Table[S] are the final photometric parameters, 
obtained by averaging our four sets of results for the individ- 
ual datasets. We found that the Xv values for these averages 
are generally noticeably greater than unity (between 1.1 and 



5 The Appendix is available only in the electronic version of this 
work. 



2.3), which suggests that there is additional uncertainty not 
picked up by our random or systematic error analyses. The 
error bars in our final photometric parameters have been in- 
creased to account for this additional uncertainty. The situ- 
ation for WASP-4 echoes that for HD 189733 and HD 209458 
l|Southworthll2008l ) and should probably be attributed to the 
effects of starspots. One corollary of this phenomenon is that 
a detailed study of a TEP cannot rest on high-precision pho- 
tometry of only a single transit: the additional uncertainty 
may then exist but not be detectable. 



4 THE PHYSICAL PROPERTIES OF WASP-4 

The parameters determined from the light curve analysis 
are not enough in themselves to calculate the physical prop- 
erties of the components of the WASP-4 system, so addi- 
tional constraints must be sought from other types of ob- 
servations and also from theoretical models of stellar evo- 
lution. The observational constraints are the velocity am- 
plitude of the parent star, Ka = 247.6^g 3 g 9 ms _1 , plus 
its spectroscopically-derived effective temperature and sur- 
face metal abund ance, T e s = 5500 ± 100 K and [^] = 
-0.03 ± 0.09 dex (|Gillon et al.l 120091 ). For t heoretical con- 



straints we use tabulations from the Claret j|Claretll2004bl . 
120051. 120061. 120071'). Y 2 l|Demarque et al.ll2004h and Padova 
l|Girardi et al.ll2000h models. 



Our approach follows that of ISouthworthl (2009): the 
above input quantities are specified and the (unknown) ve- 
locity amplitude of the planet is freely adjusted to find the 
best match between the properties of the star and the predic- 
tions of one of the sets of theoretical models. The process is 
repeated many times whilst varying each input quantity by 
its lcr uncertainty, building up a c omplete picture of the solu - 
tion and its random-error budget l|Southworth et aTll2005ah . 
There is a clear theoretical dependence arising from the use 
of stellar models in this procedure, and a lower limit on the 
resulting systematic errors are inferred by comparison be- 
tween the results for the three different sets of theoretical 
predictions. This is only a lower limit because independent 
models have some physical ingredients in common, for ex- 
ample opacity tables. As an external check we ignore the 
stellar models and instead force the star to match an empir - 
ical mass-radius relation, constructed bv lSouthworthl (|2009D 
from well-studied eclipsing binary star systems. There is one 
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Table 4. Parameters of the JKTEBOP best fits of the light curve of WASP-4 for different LD laws with the coefficients fixed at theoretically 
predicted values. For each part of the table the upper quantities are fitted parameters and the lower quantities are derived parameters. 
To is given as HJD - 2454000.0. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


rA + r b 
k 

i (deg.) 


2103+ 0021 
u.ziuo_ 0006 

1 5945+0-0004> 
U - 1O ^°-0. 00029 

89 84 +0 - 61 


2097+ 0018 

U.ZUMf _ Q 000 g 
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2097+ 0020 
— Q.0Q06 

1 A AAA^A 
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u - ldozz -0. 00025 
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1 A AAAOA 

u.iocsio_ 0()024 
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A 1 K077+O.OOOJ7 
V.IOOI 1 _ 0(H )2G 

89.73l°j 7 


"A 


0.60 fixed 


0.40 fixed 


0.25 fixed 


0.70 fixed 


0.50 fixed 


V A 




0.25 fixed 


0.50 fixed 


0.23 fixed 


0.10 fixed 


To 


36B.916894tS;ggggg 


365.916893t™47 


365.916894!°;— 


365.916894!°;— 


365.9168941°;— 


rA 
r b 


0.18251°;°™ 
0.02782_ 00011 


02790 +() """- :l 
u - uz ' MU -o. 00010 


n 1S1Q+ - 0017 
0.1819_ Q 00Q5 

0.0278b_ oooog 


o-i82o±S:SS£S 

n no7S7+ -000 35 
U - U ^' S '-0.00008 


U.U2788_ OOQOg 


a (mmag) 


0.8792 


0.8576 


0.8558 


0.8555 


0.8553 


Y 2 

A red 


1.4308 


1.3189 


1.3187 


1.3164 


1.3107 



Table 5. Parameters of the JKTEBOP best fits of the light curve of WASP-4 for different LD laws with all LD coefficients included as 
fitted parameters. For each part of the table the upper quantities are fitted parameters and the lower quantities are derived parameters. 
T is given as HJD - 2454000.0. 



Linear LD law 



Quadratic LD law Square-root LD law Logarithmic LD law 



Cubic LD law 



rA + r b 


a 9n«t;+ 0019 

u.zuoo_ 0008 


A 9ASS+ ' 0023 
U.2U8»_q 0011 


2093+ 0019 
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u.zuyi_ .o 012 


2093+ - 0022 
u.^uy^_ 0015 


k 


15471+0-00054 
U-lO<Ul_ 00041 
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i (deg.) 
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89.57+^g 
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UA 
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U.Z(0_q 2 g 4 
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0.081™ 
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u -^ lz -0.516 


114+0-isi 


093+ 144 


T 


365.916894lg : — 


365.9168851°;— 7 


365.916885_ 000050 


365.916885+2°^ 


365.916884"l°;f°°°» 



rA 
rb 

(j (mmag) 

red 



A 1SAfi+ 0015 
0.18Ub_ 0>0007 



0.1809 



+0.0019 



0.02794 



+0.00033 
0.00013 

0.8568 
1.3115 



0.02791 



0.0010 
+ 0.00053 
-0.00011 

0.8518 
1.1364 



0.1813+ Q ;^ 
0.027941SS 
0.8517 
1.1362 



A 181 9 + 0.0018 

u.iai2_ 0011 

fO. 00047 
-0.00014 

0.8517 



0.1814 



+0.0018 



0.02794^ 



0.02797 



1.1364 



0.0013 
+ 0.00053 
0.00017 

0.8518 
1.1365 



output parameter which has no dependence on the mass- 
radius relatio n or theoretical predictions: the s urface gravity 
of the planet l|Southworth et alj|2004cl . l2007bl) . 

Table[7] shows the results of the above process for 
WASP-4, taking the light curve parameters fr o m Se ction[3] 
and the K A , T e g and [M] from iGillon et al.l lj2009D . The 
variations between results for different sets of stellar models 
are most important for a and Ma, but for all parameters 
the systematics are similar to or smaller than the random 
errors. The results found using the mass-radius relation dif- 
fer by a similar amount to the systematic errors, being in 
general smaller. This probably sets an upper limit on the 
systematic errors in the properties of WASP-4, but further 
work is still necessary to illuminate the reasons why theo- 
retical mo dels persistently underestimate the radii of low- 
mass stars (|Clausenlll998l ; lTorres fc Ribadi^OOl iRibas et all 
2008). The Y 2 models yield the most precise of the age es- 
timates, due to the smooth evolution through the main se- 
quence predicted by these models, but the age of WASP-4 
remains poorly constrained. 

In Table[8] we give the final physical properties of the 
WASP-4 system, including both random and systematic er- 
ror estimates for all parameters, and compare these to those 
found by other studies. The agreement is in general very 
acceptable, even though the error bars are small. WASP-4 



therefore joins the relatively small group of best-understood 
TEPs, helped by the large amount of high-quality photo- 
metric data available for it. The error budgets we find for 
its physical properties indicate that an improved value for 
the velocity amplitude of its parent star would be useful, 
plus an improved understanding of the stellar T e s and [H] 
and the systematic errors inherent in the effective tempera- 
ture scale of low-mass stars. 



5 SUMMARY AND CONCLUSIONS 

Transiting extrasolar planets (TEPs) are vital as our pri- 
mary source of the physical properties of gas giant planets. 
The precision with which we can determine the masses and 
radii of these objects is critically dependent on the observa- 
tional data we can obtain for them, in particular the qual- 
ity of the light curves of transit events. Defocussing your 
telescope allows a much better control of systematic effects, 
many of which arise through pixel sensitivity variations, and 
also lowers the amount of observing time lost to CCD read- 
out. Telescope defocussing is rapidly becoming a standard 
method for observing TEPs, and in this work we have pre- 
sented defocussed observations of four transits of the WASP- 
4 system. 

We analysed these transit light curves in order to define 



© 0000 RAS, MNRAS 000, 000-000 



High-precision defocussed photometry of WASP-4 7 



Table 6. Final parameters of the fit to the light curves of WASP-4. The results of published studies are included for comparison, and 
are without error bars if they were not directly quoted quantities. 



Source 


This work 


This work 


This work 


This work 


This work 


Wilson et al. 


Gillon et al. 


Winn et al. 
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(Euler) 
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(Magellan) 


(20081 
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(20091 



»"A + 

k 

i (°) 
VA 



n, 
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0.00065 
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+0.0011 
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+0.00022 
0.00014 
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0.1825 
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2091+ - 0016 
"'^ -0.0025 
,+0.0012 
-0.0014 

89.6±1 A 
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u.isi2_ 0020 
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2101+ - 0053 
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2100+ - 0023 
u.ziuu_ 0022 

,+0.0012 

-0.0010 

-1+0.7 

• x -1.4 

-,+0.0021 

-0.0018 

0.02805_ Q ()m|4() 



0.154r 



0.1820 1 



0.2118 
0.15374 
88.6 

0.1836 
0.02822 



+0.0015 
-0.0017 
+0.00073 
-0.00079 
+0.6 
1.2 
+0.0012 
0.0015 
+0.00028 
0.00030 



0.2190 
u.iooz_ 00()6 
88.5911^ 
0.1895 
0.02942 



0.15353 
89.35 



0.2078 

+0.00033 
-0.00026 
+0.64 
0.49 

0.1801 



0.02764 



0.2092 

+0.00077 
-0.00055 
+0.98 
0.46 

n 1097+00017 

U.18ZV_ 00Q5 

0.027878 



0.15375 
88.56 



Table 7. Physical properties for WASP-4, derived using either an empirical stellar mass— radius relation or the predictions of different 
sets of stellar evolutionary models. Mass, radius, surface gravity and density are denoted by M, R, logg and p, respectively, subscripted 
with an 'A' for the star and a 'b' for the planet. 



Mass— radius 



Padova models 



Y 2 models 



Claret models 



a 

M A 
R A 

logSA 
PA 

M b 
flb 
Sb 
Pb 



(AU) 
(M S ) 

(Rq) 
[cgs] 

(p©) 

(M Jup ) 
(Rjup) 

(ms" 2 ) 
( Alup) 

(Gyr) 



0.022669 ±0.00058 
0.867 ±0.067 

0. 889 ±0.024 
4.478 ±0.012 

1 oqo+0.020 

1. zoo_ 022 

1 991 +0.093 
'■"-0.071 

1.334 ±0.036 



17.03 



+0.97 



0.514 



0.54 

+ 0.033 



0.022941S 
0-898±g;g!g 

900+ 020 
4 483+0 011 

1 9 qn + 0.020 

i.zoo_ 022 

1 ncrn + O.OUO 

i.zou_ 056 

1 ocn + 0.032 

i.oou_ 025 

17.03^-JI 

0.508±g : °il 
7 0+2.3 

'■ -4.2 



n no'RQn+O-OO 036 
u.uzjzu_ () 0()()4H 

009+0 032 

910+ 015 
u.»iu_ 020 

4 488+ ' 008 

1 9 qn + 0.020 

i.zoj_ Q22 

1 O7Q+0.082 
l.ZI3_ 063 

1 ocK+0.023 
±.ouj_ 02g 

17 Q3+0.97 
i(.uo_ 54 

502+ 031 
u.ouz_ 02Q 

5 5+ 3 ' 2 
°'°-2.0 



02329+ - 00040 
Q40+0 048 
u,3 * u -0 .069 

014+O OI6 
u.yit_ 023 

4 490+ 009 

'4.1»U_ Q12 

1 9 oq+0.020 
i.zoo_ 022 

1 9SQ +0.085 
l.zc.y_ Q73 

1 071+O.O26 

J^-^^-o.oss 

17 Q3+0.97 
i(.uo_ 54 

500+ 032 
u.ouu_ 020 

R 9 +53 



the photometric parameters for WASP-4 (ta, r^, i and pe- 
riod) , with careful attention paid to the incorporation of stel- 
lar limb darkening and to the assessment of random and sys- 
tematic errors. Three other sets of high-quality light curves 
(precision better than 2 mmag per observation) exist in the 
litera ture (jWilson et al.ll2008l : lGillon et al.ll2009l : IWinn et~al1 
2009), which we have modelled in the same way. Combining 
these analyses with our own maximises the quality of the fi- 
nal results and minimises the possibility of systematic errors 
caused by starspots or instrumental effects. The photometric 
parameters were then added to the observed stellar veloc- 
ity amplitude, effective temperature and metal abundance 
in order to determine the physical properties of the star and 
planet in the WASP-4 system. We undertook this step using 
additional constraints from theoretical stellar evolutionary 
models, and assessed the resulting systematic error by com- 
paring results obtained using different sets of models. 

We find that WASP-4 b has 

Mj up , 



1 OSO. +0 090 +0.039 
l.ZOa -0.073 -0.000 



a mass of 
radius of 



R h = 1.371 "tn°n 3 -! 2 s Innnn Rjup, and a surface gravity 



of gb 



17.03 



0.035 -0.000 
+0.97 „„-2 



ms .Its equilibrium temperature 



with the trend s for short-pe r iod p lanets to have relatively 
large masses ( Mazeh et al. 2005) and surface gravities 
|Southworth et aL 2007bl ). The oretical models of coreless 
irradiated gas giant planets (p3odenheimer et all 120031; 
iFortnev et all 120071 ; iBaraffe et aL 200Sl ) give predictions 
of the radius of WASP-4 b ranging from roughly 1.09 to 
1.17 Rjup, and are therefore rejected at beyond the 5a 
level. Models with a heavy-element core propose smaller 
radii and are therefore even more discrepant. However, 



alternative mo dels including an ad hoc additional kinetic 
heating source l|Bodenheimer et alJliooll ) can satisfactorily 
explain the properties of WASP-4 b. 

The high quality of the available data for WASP-4, 
plus the comparative depth of the transits (3%), means 
the mass and radius measurements of its planet are 
among the most precise a vailable, toge ther with C0R0T - 
Exo-2 (lAlonso et alll2008l). H P 18973 3 JPont et al.ll2007l) . 
HD 209458 (ISouthworthll2009T) TrE S-3 (|Sozzetti et alj|2009h 
and WASP- 10 (| Johnson et al.ll2009l ). Error budgets indicate 
that a more precise velocity amplitude for the parent star 
would be useful for better measuring the physical properties 
of WASP-4 b, but further improvement will require tackling 
the systematic errors which are unavoidable in the analysis 
of TEPs. The predictions of different sets of theoretical stel- 
lar models simply disagree for low mass stars, both between 
themselves and with the observed properties of low-mass 
eclipsing binary star systems. In the case of WASP-4, this 
disagreement is at the 5% level for the mass of the planet 
and the 8% level for the star's mass, with other properties 
less affected. This level of disagreement is a problem for in- 
ferring the chemical com positions of gas gian t planets from 
their masses and radii (see lFortnev et aLll2007t ). For the best- 
studied cases, our understanding of TEPs therefore remains 
limited by our lack of understanding of low-mass stars. 
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Tabl e 8. Final physical p roperties for WASP -4, B oth statistical and systematic uncertainties are included. The corresponding results 
from IWilson et all l l2008t) , iGillon et ail feOOSft and IWinn et al.l j2009l) have been inclu ded for comparison. The surfac e gravity of the 
planet has no systematic error as it can be specified by purely observational quantities l|Southworth et al , 2004c, 20 7bh . 



Final result (this work) Wilson et al. (20081 Gillon et al. (2009) Winn et al. (2009) 
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Ra 
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Rh 
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Pb 



(AU) 
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(ms" 2 ) 

(Alup) 



n rwwri +0.00050 +0.00058 
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Table Al. Parameters of the JKTEBOP best fits of the Euler i?-band light curve of WASP-4 [|Wilson et al.l2008h . using different approaches 
to LD. For each part of the table the upper quantities are fitted parameters and the lower quantities are derived parameters. To is given 
as HJD - 2454000.0. The light curve contains 213 datapoints. 
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n Q81+0 034 
u.ool_ 038 


VA 




0.25 perturbed 


0.50 perturbed 


0.23 perturbed 


0.10 perturbed 


To 




365.9159791°;— 


365 Q15Q75+0-000 110 
' 3DO - yl09 ' -0. 000103 


365 015077+0.000101 
■soo.yioyr <_ . 000104 


365.9159751°;™ 


VA 
rb 




n TS19 + 0.00 32 

0.1812_ Q014 
0.0285b_ 00021 


n 1815+0 00 36 
0.1815_ 0016 

n n9a79+o |,,,,,7i; 

O.U2»|-2_ 00026 


n 1814+ - 0034 
0.1814_ 0Qig 

n n98fi5+0- 00078 

U - U/8bt »-0.00024 


0.18151™ 
n n9S75+ -00 082 

U.02SC5_ 00030 


a (mmag) 




1.9261 


1.9254 


1.9255 


1.9253 


X red 




3.2553 


3.2527 


3.2533 


3.2523 


Fitting for both LD coefficients 


rA + r h 
k 

i (deg.) 
u A 

VA 

To 




2100+ 0046 
u.ziuu_ 0025 

158fi2+ ()l)017 ' s 
U - 10SDZ -0. 00163 

88.8+J- 3 

0.370lg;'^? 

0.069 + o 2 53 
365 01507+0.00011 

000.910M(_ Q 00012 


2102+ - 0040 
u.ziuz_ 0029 

^58fi0+ () ■ l)0 '- i " , 
U.lOSOU_ 00227 

88.8+4-5 

252+ - 715 

259+ 1 - 243 
u.zoy_j 243 

365 01507+0.00011 

■ 5D0 - a101 " -0.00011 


2101+ - 0039 
u.ziui_ 0026 

15855+0.00196 
U.lOOOO_ 00191 

88.8+4-5 

47 l+0.'269 

u -*' -0.267 

097+ - 336 
u - ua ' -0.358 

365 91597+0 000 11 
dDO.9ioy/_ 00011 


2102+ 0049 
U.Z1UZ_ 0023 

fl 1 5840+ ' 1 " 12112 
U.1584y_ Q 00197 

88.8lj-g 

n Q89 + °'0 66 
u.ooz_ 06g 

090+ 0,281 
U.U3U_ 295 

365 91597+0 00011 
csDO.9ioyi-_ 00011 


rA 
rb 




n ioiQ+0.0036 
u - lolo -0 .0022 

n n9«7fi+o ° (lo: ' :i 
0.0287b_ 00031 


fl 1815+0-00 32 
U.±O±O_ 0027 

n n9S7a+o ooo')."i 

U - U ^ 8 ' 8 -0.00029 


n 1S14+ - 0033 
0.1814_ 0022 

n n9a75+0' ()0H91 

U.UZcWO_ Q 00030 


n 181 5+0.0039 
u.ioio_ 0020 

n n987fi+ H n 

0.0287b_ 00027 


a (mmag) 




1.9253 


1.9254 


1.9253 


1.9253 


^ red 




3.2682 


3.2683 


3.2682 


3.2681 



APPENDIX A: RESULTS OF THE LIGHT 
CURVE ANALYSES 

The tables in this section contain the full results of modelling 
each light curve of WASP-4. 
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Table A2. Parameters of the JKTEBOP best fits of the z-band VLT light curve of WASP-4 l|Gillon et alJ l2009h. using different approaches 
to LD. For each part of the table the upper quantities are fitted parameters and the lower quantities are derived parameters. To is given 
as HJD - 2454000.0. The light curves contain 244 datapoints. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


All LD coefficients fixed 


rA + r b 
k 

i (deg.) 


208 c i +a0017 
u.zuoo_ 0004 

U.lOlOI_ q 0()23 

SQ 01+0 46 
oa.yi_o g5 


2099+ 0015 
u.zuyy_ 0014 

0.15304_°;°°™ 

88.561°;™ 


2094+ 0016 
u.zuy4_ 0015 

(1 1 cqi o+0. 00027 
u.uoio_ .00026 

88.61_"°;^ 


2094+ 0015 
u.zuy*_ 0015 

0.153081HS 
88.64±°;£ 


2092+ - 0016 
u.zuyz_ 0015 

1KQ1Q+0-00031 
U.lOOlO_0 oo()2« 

88.65.0^ 


ua 


0.50 fixed 


0.25 fixed 


0.10 fixed 


0.59 fixed 


0.40 fixed 


VA 




0.31 fixed 


0.54 fixed 


0.26 fixed 


0.10 fixed 


To 


396 695401+ 000037 

OMD.D»O1Ul_ .000036 


396 695391+ 000035 
oau.u»uo»±_ .000035 


396 695394+ 000036 
oau.uajoi>»_Q 000035 


396 695393+ 000036 
o»u.uyooyo_o 000033 


396 695395+ 000035 
oyu.oyooyo_ 000037 


rA 
rb 


n i8in+ a0014 
u.ieiu_o.ooo4 

027451+ - 000271 

U.UZI401_ Q 000076 


n 1090+0.0013 

0.1820_ 0012 

n n97Rc:+0' l)002:i 
U.U2C8b_ 00023 


0.1816_°;°£ 3 3 
0.02780lg;°r2 6 3 


0.181b_ 0013 
fl p,97S1+0'00024 
U.U2781_ 00023 


n 1S14+0- 0013 

0.1814_ 0013 

n n977«+o | »i»)2^ 

UM2/ib_ 00024 


a (mmag) 


0.6650 


0.6111 


0.6056 


0.6059 


0.6075 


Y 2 

^ red 


2.0671 


1.6927 


1.6552 


1.6571 


1.6680 


Fitting for the linear LD coefficient and fixing the nonlinear LD coefficient 


TA + r b 
k 

i (deg.) 

UA 


2103+ 0016 
U.ZiUO_o 0016 

0.15442«S 
87.96+*f 9 

Q ogn + 0.010 

u.oyz_ Q 011 


0.208b_o 0007 
15244+"-"" 014 

U.ldZ^_ 00033 

89 42+ - 66 

270+ 013 
u - z,u -0.013 


2099+ 0015 
u.zu»y_ 0017 
fl icocri +0.00042 

00 qqHtO.63 
88.33_o 47 


2094+ ' 0017 
u.zuy*_ 0014 

0-15307_°;°°°« 

88.65_°;g 

590+ 012 
u.oyu_ 011 


2099+ - 0016 
u.zuyj_ Q 001s 

0.15378l 'Z 44 

88.25l0- 6 44 

n o 7 c:+0.bl2 
U.JIO_o 012 


VA 




0.31 fixed 


0.54 fixed 


0.26 fixed 


0.10 fixed 


To 


396 695394+ 000034 
o»u.u»jo»i_ .000035 


396 695393+ 000035 
oau.u»uo»o_ . 000034 


396 695393+ 000038 
oyD.oyooyo_ 000036 


396 695393+ 000034 
oyu.uyjoyo_ 000037 


396 695393+ 000035 
jyD.oyojyo__ . 000036 


rA 
rb 


n 1091 +o.ooi3 

U.18^!l_o 0013 

n n9si-?+o °oo 26 

U.lUcUd_ 00026 


n i8in+ ' 0014 

U.l»10_o 0006 

02759+ - 00028 
U.UZfO»_ 00010 


0.1820_ 0014 

02793+ 00027 
u.uzi ao_ .00027 


n isin+ a0014 
0.181b_ 0012 

n n97sn+ lMlll(l2<l 
0.02780_ 00023 


n isiq+ - 0014 

0.1819_ 0015 

02797+ - 00027 
u.uziyi _ Q0029 


a (mmag) 


0.6046 


0.6095 


0.6049 


0.6059 


0.6046 


red 


1.6533 


1.6900 


1.6566 


1.6640 


1.6538 


Fitting for the linear LD coefficient and perturbing the nonlinear LD coefficient 


ta + r b 
k 

i (deg.) 
u A 




U.^08b_ .ooo8 

15244+ - 00055 

U.lOZ^_ 00048 

89 42+0-71 

oy.4z_o.8 4 

270+ ' 027 
U.Z(U_ 029 


2099+ 0016 
u.zuyy_ 0016 

fl 1W1 +0.00045 
U. IOOOI _ oqo.39 
00 qq+0.65 
88.33_o 47 

o.oselHfg 


2094+ 0018 
u.zuy*_ 0013 

1 5^n7+ - 00057 

U.lOOU(_ .00049 

590+ 051 
u.oyu_ 053 


2099+ 0015 
u.zuyy_ 0017 
n 1 KQ 7 o+0. 0001s 
U.!OOIO_o .00000 

88.25ir4 8 5 

n o 7 c:+0.016 

u.jio_ 016 


VA 




0.31 perturbed 


0.54 perturbed 


0.26 perturbed 


0.10 perturbed 


To 




396 695393+ 000037 
o»o.uyoojo_ Q .000036 


396 695393+ 000033 
oyu . oy ooyo _ .000034 


396 695393+ 000038 
oyu.oyooyo_ 000032 


396 695393+ 000034 
oyo.uyooyo_ 000037 


rA 
rb 




0.1810_o 0007 

02759 +(K0 "" :!0 
u.uzioy_ 00010 


n 1 o9n+0'0 013 

0.1820_o 0014 

02793+ 04 27 
u.uzi ao_ 00028 


0.181b_o 0011 

n n97sn+ [, ■ < " lll2!, 
0.02780_ 00024 


n isiq+O- 0013 

0.1819_ 0014 
n 619707+° .0UU26 
u.uziyi _ 00030 


a (mmag) 




0.6095 


0.6049 


0.6059 


0.6046 


Y 2 

red 




1.6900 


1.6566 


1.6640 


1.6538 



Fitting for both LD coefficients 



rA + r b 
k 

i (deg.) 

UA 

VA 

To 


2100+ 0017 
U.ZlUU_o 0017 

15412+ 000067 
U.lO41Z_ 00076 

88.12l°; 72 

n q71 +0.037 

u - J,± -0.O39 

050+ - 091 
U.UOU_ o 8 3 

396 695393+ 000034 
oyo.uyooyo_ 000035 


2101+ 0016 
U.Z1U1_ 0016 

15412+ 041 ""*" 
u.lO'ilz_o 00089 

88.09l°-™ 

290+ ' 250 
U.ZJU_ 248 

179+O.425 

u,lla -0.440 

396 695393+ 000033 
oyu.uyooyo_o 000035 


2100+ 0017 
U.Z1UU_ 0015 

15412+ - 00078 

U.lOtlz_ 00()7 3 

00 11+0.69 
88.ll_o.4o 

Q 4 o ? +0.091 
"■^1 -0.097 

0.059 + oi23 

396 695393+' 000036 
oyu.oyooyo_ 000034 


2101+ - 0017 
U.Z1U1_ 0017 

88.08l°; 74 

n q«t;+0.023 
U.OOO_ Q23 

n ,: i9+ 0115 
u.uoy_ 099 

396 695393+ 000036 

030.030030_ 000038 


rA 
rb 

a (mmag) 

Y 2 

^ red 


0.1819_ 0014 

n n9snzL+0' H00:i2 
0.02804_ 00032 

0.6042 
1.6586 


U.18Z^U_o 0013 
n n98rK+ a °oo3i 

O.UZ^80o_o.oo035 

0.6044 
1.6596 


U.18^U_ .ooi3 

02805+ ' <lllll:u 
u.uzeuo_o .00030 

0.6044 

1.6592 


0.1820±g;ggli 
0.02806l°;^ 32 
0.6044 
1.6595 
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Table A3. Parameters of the JKTEBOP best fits of the Magellan z-band telescope light curve of WASP-4 l|Winn et al J 12009V using 
different approaches to LD. For each part of the table the upper quantities are fitted parameters and the lower quantities are derived 
parameters. To is given as HJD — 2454000.0. The light curves contain 713 datapoints. 



Linear LD law Quadratic LD law Square-root LD law Logarithmic LD law Cubic LD law 



All LD coefficients fixed 



k 

i (deg.) 
ua 
va 
To 



+0.0014 
0.0004 
+0.00028 
-0.00021 
+0.56 
0.79 
0.50 fixed 



0.2103 
0.15299 
89.63 



697 7O7572+0-000032 
DM(.(»IO/Z_ 000038 



0.2126: 



,+0.0014 
3 -0.0014 

U.104UO_ 00025 

88.26l j 7 
0.25 fixed 
0.31 fixed 



697.797565 



+0.000034 



0.2122 
0.15428 



+0.0015 
0.0013 
+0.00030 
-0.00025 
+0.45 



88-25_ .4o 
0.10 fixed 
0.54 fixed 
697 797567+0-000032 
— 0.000036 



2124+ 0014 
U.ZlZ4_o 0016 
15419+0 00026 
U.lJ<iiy_o .00029 

88.24l°jl 
0.59 fixed 
0.26 fixed 
697 797566+ 000033 

DM(.(»/ODD_ 000035 



+0.0015 
0.0013 
+0.00028 
-0.00030 
+0.50 
0.40 
0.40 fixed 

0.10 fixed 



0.2119 
0.15431 
88.31 



697.797568 



+0.000036 



n ,09^+0.0012 n 18ZL9+ ' 0012 n 1 sqq+o 0012 ni84n+° oo 12 n i83fi+ - 0012 

r A 0.1824_ooo04 U - 1842 -0.0012 °- 18iJ9 -0.0011 U - 184U -0.0013 °- 18db -0.0Qll 

n n97Qn+°-000 24 n r^s^a+O- 00023 n n9s-?7+o ooo 24 n rws^s+O- 00021 n rwR^+O- 00023 

r b O.U2|-yO_ 00008 U - U28d8 -0.00022 U.lWBiU _ .00021 U - U28 ' ;i8 -0.00024 U - U28 ' ;i ' ;i -0.00022 

cr (mmag) 0.6479 0.6039 0.6013 0.6012 0.6027 

X 2 rod 0.7892 0.7003 0.6944 0.6947 0.6964 



Fitting for the linear LD coefficient and fixing the nonlinear LD coefficient 



rA + r h 
k 

i (deg.) 
ua 
va 
To 



0.2123 



-0.0015 



-0.0014 

o i r )->o( : r""" 0:! ' ! 

U.100UO_Q 0003g 
-,+0.011 



0.412^ 



0.012 



697.797568 1 



,+0.0016 
-0.0009 
15S^ fi +0 00046 

88.98*" 
0.276+°°" 



0.2112 1 



697.797565 



0.014 
0.31 fixed 

+0.000036 
0.000033 



-,+0.0016 
-0.0016 
,+0.00041 
-0.00046 

88.27+°" 

,+0.012 
-0.013 
0.54 fixed 

697 797567+0 00 0034 
DM '-' y ' 00 '-0.000034 



0.2122Z 
0.15424 4 



0.102 4 



3+0.0015 
5 -0.0014 
fl 1 cooy+0. 00044 
U - 10,30 '-0.00041 
,+0.70 
-0.50 
5+0.012 
-0.013 

0.26 fixed 
697 797566+°-° 00034 

DM ( . (»IODD_o .000032 



0.2118: 
15387 4 
88.51 H 
0.603 4 



n 9199+0.0014 

u.zizz_ 0015 

1 5452+0-000 42 
U.1040^_ 00041 

88.18l°;f 
391+0.011 
u.oa±_ 012 

0.10 fixed 

+0.000034 



697.797567 



-0.000034 



rA 1838+ ' 0012 
r A U - i8 ' ; ' 8 -0.0012 

02849+ - 00025 
r b u.uzo4»_ 00023 

<7 (mmag) 0.6039 

X 2 rcd 0.7014 


n 18^1 +0.0013 
u. 100 -0.0Q07 

0.02808+0 1 ™ 
0.6020 
0.6961 


n ioqo+0. 0013 
0.1S38_o 0013 

0.028361°;°°°!? 
0.6014 
0.6954 


n 1S o fi +o.ooi3 

0.1836_ 0012 

02825+°-°°° 26 
U.UZ8ZO_o .00024 

0.6009 
0.6943 


n ioqo+0.0012 
0.1»3»_ 0013 

n ri98 in+o-ooo2.-, 
0.02840_ Q 00025 

0.6018 

0.6965 


Fitting 


for the linear LD coefficient and perturbing the nonlinear LD coefficient 






rA + r b 
k 

i (deg.) 
u A 

VA 

To 




n 911 9+O.OOI6 
U.Z11Z_ 0009 

n i coo fi +o.ooo60 

U.lOOOD_ 00040 

88.98+°-^ 

276+ 030 
U.Z(D_ 029 

0.31 perturbed 

697 797565+ 000033 
mt .<viovo_ ooo36 


n 9199+0.0015 

u.zizz_o.ooi6 

15424+0-0 0047 
U.lO^Z4_ 00040 

88.27+°- 4 2 
102+ 038 

U.lUZ_o 039 

0.54 perturbed 

697 797567+0 0000 33 
— 0.000032 


n 91 10+0. 0016 

U.211S_o 0015 
fl 1 eqo 7 +0. 00053 
U.iOOO/_o .00054 

88.511°;^ 

0.603lg;gg 

0.26 perturbed 

697 797566+ 000033 
' • ' " ' — 0.000033 


n 9199+0.0015 

U.Z1ZZ_ 0016 
0-154521°;°°°^ 

88.18ir 4 2 3 
ogi+0.016 
u.oa±_ 018 

0.10 perturbed 

697 797567+0 0000 32 

ot)(.l»IOO(_ ooo34 



n 1001 +0.0014 n ,ooo+0.0012 ni83fi+°'00 14 O 

r A 0.1S31_oooo8 Ua8d8 -0.Q013 U - i8do -Q.0Q12 U - 18d8 -0.0014 

n rwsns+O- 00031 n rws.fi+o 00027 n n989=;+o ooo 27 n n98An+0- 00025 

r b U - U28U8 -0.00015 U ' U28 ' ;iD -0.00026 U - U282o -0.00029 U - U284U -0.00029 

cr (mmag) 0.6020 0.6014 0.6009 0.6018 

X 2 rod 0.6961 0.6954 0.6943 0.6965 



Fitting for both LD coefficients 



'b 



rA 
k 

i (deg.) 

"A 
Va 
To 



2116+ 0016 
U.Z11D_ 0014 

15400+°' 00071 

u.±d<iuu_ 00066 
CQ C9 + 1.03 

32+°' 04 
u " 3z -0.04 

202+ 098 
U.zuz_ 094 

697 797566+ 000034 

' • ' a ' ow — 0.000034 



0.2119 
0.15356. 
88.60 
-0.12 
0.928 
697.797566 



+0.0017 
0.0012 
+0.00090 
0.00086 
-1.23 
-0.66 
+0.25 
0.26 
+0.457 
0.438 
+0.000035 
0.000034 



n 91 1 o+u.uu-w 
U.2118_o 0013 

fl 1 KQ7R + 0.0008!) 
U.lDofO _ . 00077 



0.2118 



+0.0016 



7+1.16 
-0.65 
,+0.10 
-0.10 
,+0.137 
-0.134 

697 797566+ 000033 
ua; . i»iouv_ Q00034 



88.57Z 
0.62j 
0.282 4 



-0.0013 
,+0.00093 
-0.00086 
hi. 16 
-0.65 
7+0.02 
-0.02 
240+ 130 

697 797566+ 000035 

o»l./a<dOU„ 000036 



0.15363: 
88.6i: 
0.37 4 " 



n 1001+O.OOI3 n ioo 7 +0.0014 +0.0014 n ioo fi +0.0012 

r A 0.1834_ooon O-i^'-O.OOlO u - 18 - J5 -o.0011 0.183b_ 00 oii 

n n9894+0-000 31 n 09091 +0.00035 n 09899+0' 00033 n 09091 +0.00032 

r b 0.02824_ 00026 °- 02821 -0.00025 °- 02822 -0.00027 °- 02821 -0.00026 

cr (mmag) 0.6009 0.6010 0.6009 0.6010 

X 2 rod 0.6952 0.6953 0.6952 0.6952 
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Table A4. Parameters of the JKTEBOP best fits of the Danish telescope ij-band light curve of WASP-4 (this work), using different 
approaches to LD. For each part of the table the upper quantities are fitted parameters and the lower quantities are derived parameters. 
To is given as HJD — 2454000.0. The light curves contain 452 datapoints. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


All LD coefficients fixed 


ta + r b 
k 

i (deg.) 

"A 


2103+ 0021 
u.ziuo_ () 0006 

15245+ ' 00045 
U.lOZto_ . 00029 

89.84+0 ^ 
0.60 fixed 


2097+ 0018 
u - zuy ' -0.0006 

15345+OOOO 33 
U.lOO1O_ 000 23 

89.82+. °j? 
0.40 fixed 


2097+ ' 0020 
u.zuy/_ oo6 

U.100ZZ_o. 00025 
SO 07+0-52 

oy.y 1 _ 1 04 
0.25 fixed 


2099+ ' 0021 
u.zuyy_ .ooo5 
1531R+000039 

U.10010_o. 00024 

89.971?; 4 ? 
0.70 fixed 


2092+°°° 2 ° 
U.ZU»Z_ 0005 

fl 1 5Q 77 +0.00037 
U.10CW 1 _ Q 00026 

89.73+°;^ 
0.50 fixed 


VA 

To 


365.916894lg;gg~g 


0.25 fixed 
365.916893lS-.gg~S 


0.50 fixed 
365.916894 4 :° ; S 44 


0.23 fixed 
365.916894 4 :° ; ~ 


0.10 fixed 

365 Q16894+ 000043 
ODO.yiDe>y^_ .000046 


TA 

n> 

a (mmag) 


0.1825lHSo 8 5 
0.02782_ ooon 
0.8792 


0-i8l8±°;gg« 

02790+ 00029 

U.UZ(MU_ 00010 

0.8576 


0.02786+ 9 ;^ 3 09 
0.8558 


0.18201°;°^ 
n n<)7S7+0'00» : ',5 
U-U2C8C _ .00008 

0.8555 


n 1R1Q+0.0016 

U.l»13_ 0004 
09700+0. 00032 

0.8553 


V 2 

A red 


1.4308 


1.3189 


1.3187 


1.3164 


1.3107 


Fitting for the linear LD coefficient and fixing the nonlinear LD coefficient 


ta + r b 
k 

i (deg.) 

MA 
VA 

To 


2085+ 0019 

U.ZUOO_ ooos 

0-15471+°;^ 
89.30+°;™ 

u.oui_ 013 
365 916894+ 000039 

ODO.yiDOMt_ 000047 


2099+ 0021 
X5344+0-00066 
u.loo**_ .00045 

89.60+°/^ 
401+ 0014 
0.25 fixed 

365.916893+™44 


2093+ 0020 
n 15009+0.00052 

U.lOO»z_ 00045 

89.50l°; 72 

220+ ' 014 

0.50 fixed 

365 Q16894+ ' 000044 
ODO.yiooyt_ oooo46 


2094+ ' 0020 

u,iTO, -0 .0007 
15363+ 00054 

89.95 + ?-®| 
0.6781H1 4 
0.23 fixed 

365.916893+.°,;™ 


n 9nss+°-00l 9 
0.2U88_ 0007 

i5407+ -00 058 

U.1J*U( _ 00041 

90.00+°'j2 
0.485_ 015 
0.10 fixed 
365.916894+™ 47 


rA 
rb 

a (mmag) 


n i«nfi+0- 0015 
0.180b_ 00 oo7 

o 02794+0000 33 

u.u^ 1 »<±_ .00013 

0.8568 


1819+ 0017 

U.101J_„ 0007 

02792+ ' 00041 
u.u^l»^_ 00010 

0.8575 


U.1814_ oo7 

02791 +0 °0 037 
u.uziyi_ . oooil 

0.8559 


n isi 5+0-0016 

U.±010_o 0007 

n n97s«+ - 00037 
o.U2/»a_ oooio 

0.8562 


0.18091°;°^ 
n no7«7+0-00036 
U.U2C»C_ 000O8 

0.8559 


V 2 

A red 


1.3115 


1.3222 


1.3116 


1.3138 


1.3109 


Fitting for the linear LD coefficient and perturbing the nonlinear LD coefficient 


ta + r h 
k 

i (deg.) 

«A 
VA 

To 




2099+ 0012 
u.^u»y_ 0009 

0.15344^0.00064 

89.60+° f s 

401+ 029 
u -* ul -0 .028 

0.25 perturbed 

ock mcooo+0. 000045 

365.916893+ () 000046 


2093+ 0020 
u.zuyo_ ooos 

15302+ '" ,I1I( ' 1 
u.iojyz_ oo43 

89.50l°;» 

220+ 038 

0.50 perturbed 

365 916894+ 000046 
JOO.MlOO»4 : _ 000043 


2094+ - 0020 
u.^uyi_ .0008 

15363+ 00064 
U.iOODO_o 00052 

89 q5+°- 58 
oy.yo_j oo 

0.678i°;°g 
0.23 perturbed 
ock 01 cono+0. 000047 

365.916893^0.000046 


n 9nss+ -°oi 8 
0-^088-0.0008 

i5407+°-0 007 

U.1041U( _ 00058 

90.00 + °f 2 

0.485_ oi 8 
0.10 perturbed 
365 916894+ 000042 

000. 310031_ .000045 


rA 
rb 

a (mmag) 




1819+ 00016 
n 09799+ 000040 

U.UZ(MZ_ 00011 

0.8575 


fl 1R14+0 00 16 
0.1814_ Q Q007 

09701+0.00037 
u.u^i yi_ oooiO 

0.8559 


1815+ 0016 
U.1010_o oo8 

n n97ss+°- (K)037 

U.U2C8»_ ooo9 

0.8562 


1809+ 00015 
u. iouy_ ooos 

n no7«7+°-000 36 
0-02 CSC _ . 00009 

0.8559 


V 2 

A red 




1.3222 


1.3116 


1.3138 


1.3109 


Fitting for both LD coefficients 


rA + r b 
k 

i (deg.) 

"A 

To 




0.2088±°;°°3 

0.15428lH°oo°7 
90.00 + ?'|° 

476+ 046 
U.*(D_ 047 

0.081+°,^ 
365.916885+°;— 


90Q3+ 0019 
U.ZUy3_ 0014 

15409+ ' 00116 
u.id^uy_ 00103 

89.58+°'^ 

276+ ' 291 
U.ZID_o 294 

4i2 +0 ' 531 

365.916885+S°°« 


n 9001+0.0022 
u.zuyi_ 0012 

15423+ 00123 
u - lo * zo -0. 00086 

89.57+?-^ 
505+0.117 

H4+0-151 
365.916885+°f°°°^ 


2093+ 0022 
u.zuyo_ 0015 

15426+ 00130 
U.lO4ZD_ 00106 

89.36t°- 92 
493+0.028 
u.tyo_ 028 

0.093+ 91 , 44 
365.916884+°;— 


rA 
rb 

a (mmag) 




0.18091°;°°}° 

02791+" 
u.uz<y±_ 00011 

0.8518 


0.1813t°;°™ 

02794+ ' 00047 
u.uz/y^_ oooil 

0.8517 


0.1812±g : ggl? 

n97Q4+ - 00047 
O.U2794_ Q ooi4 

0.8517 


0.1814+°;°^ 

097Q7+ -000 53 
u - uz,a ' -0.00017 

0.8518 


V 2 

A red 




1.1364 


1.1362 


1.1364 


1.1365 
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